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Motivated by recent ultrafast pump-probe experiments on high-temperature superconductors, we 
discuss the transient dynamics of a d- wave BCS model after a quantum quench of the interaction 
parameter. We find that the existence of gap nodes, with the associated nodal quasiparticles, 
introduces a decay channel which makes the dynamics much faster than in the conventional s-wave 
model. For every value of the quench parameter, the superconducting gap rapidly converges to a 
stationary value smaller than the one at equilibrium. Using a sudden approximation for the gap 
dynamics, we find an analytical expression for the reduction of spectral weight close to the nodes, 
which is in qualitative agreement with recent experiments. 
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Introduction. - Recent advances in time-resolved spec¬ 
troscopies have triggered a growing interest in the tran¬ 
sient dynamical behavior of high-temperature supercon¬ 
ductors optically excited far from equilibrium. By shin¬ 
ing the sample with intense ultrafast pulses (pump) one 
can trigger nonequilibrium transient states, whose physi¬ 
cal properties are then recorded by a second pulse (probe) 
which hits the sample with a given time delay. This ap¬ 
proach opens up a wealth of information unavailable to 
conventional time-averaged spectroscopies [1-9]. When 
used in combination with angle-resolved photoemission 
spectroscopy (ARPES), these ultrafast methods allow us 
to track and follow in real time the evolution of quasipar¬ 
ticle modes in different momentum sectors, such as those 
close to the gap nodes and antinodes of d -wave cuprate 
superconductors [10-12]. In addition, when irradiation 
is sufficiently strong and optically tailored in such a way 
to selectively excite specific modes, one can even stabi¬ 
lize transient states with fundamentally different physi¬ 
cal properties [13-15]. A spectacular example is given by 
recent experiments reporting signatures of light-induced 
metastable superconductivity in cuprates at much higher 
temperature than at equilibrium [16, 17]. 

These experimental breakthroughs raise a number of 
intriguing questions. From one side, a nontrivial and 
rich transient dynamical behavior is expected in corre¬ 
lated materials, which feature complex phase diagrams 
characterized by competing phases [18]. On a more 
fundamental level, pump-probe experiments suggest the 
possibility to explore novel metastable phases that can 
only be accessed along nonthermal pathways, e.g., by 
means of photoexcitation. Theoretical investigations 
along these directions have appeared in the literature 
in recent years [19-22] addressing questions like the 
thermalization of pump-excited Mott or Kondo insula¬ 
tors [23, 24] and the role of lattice vibrations [25], orbital 
degrees of freedom [26] and competing orders [27, 28] in 
the relaxation dynamics. The dynamics of conventional 
superconductors has also been studied with reference to 


pump-probe experiments [29-31] or in the presence of 
electron-phonon interactions [32-34] and Coulomb repul¬ 
sion [35]. In the context of ultracold atomic systems, the 
research has focused on the s-wave weak coupling BCS 
regime [36-39], with recent attempts to extend the anal¬ 
ysis to the crossover into the BEC regime [40, 41] and to 
exotic order-parameter symmetries [42, 43]. 

Yet, despite this recent activity, many fundamental 
questions concerning specifically the transient response 
of superconducting materials remain wide open. In par¬ 
ticular, a characteristic feature of high-T c cuprates is the 
momentum anisotropy of the superconducting gap, which 
has a d x 2_ y 2 symmetry (which we indicate as d -wave in 
the following). This leads to the existence of nodal lines 
along which the superconducting gap vanishes, a charac¬ 
teristic feature which is sometimes used as a definition of 
exotic superconductivity. These gapless excitations sig¬ 
nificantly affect the thermodynamics and spectral prop¬ 
erties as compared to the conventional s-wave case, where 
the gap is uniform over the whole Fermi surface, and they 
can possibly play an even bigger role in the nonequilib¬ 
rium dynamics. 

In this Letter we consider a minimal model of a d -wave 
superconductor and the simplest nonequilibrium proto¬ 
col of a quantum quench of the interaction parameter. 
Within a BCS-like mean-field approximation, we calcu¬ 
late the gap dynamics and compare the results with those 
for the conventional s-wave superconductor, identifying 
the distinctive features descending from the existence of 
the nodes. Then, we derive an approximate formula for 
the spectral weight which is directly relevant to ARPES 
experiments in d -wave superconductors. Of course a 
mean-field description, while reasonable for the deeply 
overdoped region of the phase diagram of cuprates, is un¬ 
able to capture the strong correlations which dominate 
the low-doping region, leading to a complex interplay be¬ 
tween different phases. However, we believe it is useful to 
highlight the physics of a simple model with d -wave sym¬ 
metry, as this will help to disentangle the effect of the 
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anisotropic gap in more involved calculations in which 
competing orders, strong-coupling effects, and any other 
realistic feature are included. 

Model and gap dynamics. - As a starting point of our 
analysis we consider the two-dimensional BCS Hamilto¬ 
nian with a momentum-dependent separable interaction 
in the d- wave channel 

H = c L 7k7 P 4 t c -i 4 c _pi c pt , (1) 

ktr k p 

with e k = |k | 2 — p, and y k = cos 2 6 ~ k 2 — k 2 where 8 is 
the polar angle in k space. The superconducting gap, or 
order parameter, is defined as 

A k = A7 k , A = J ^2 7 k (c-k4.c k f), (2) 

k 

and it vanishes along the nodal lines k x = ztk y . 
This leads to modes at arbitrary low energies E k = 
y/e k + |A“ | which, as we show below, have important 
consequences on the relaxation dynamics. 

The existence of the nodal lines marks a fundamental 
difference not only with the s-wave symmetry which has 
7 k = 1 and a fully gapped spectrum [37, 39], but also with 
unconventional symmetries such as p + ip ( 7 k ~ e“) [42] 
and d+id (y k ~ e 2 ™) [44] where the gap vanishes at most 
for one point in momentum space. Notably, in all the 
above nodeless cases the 9 dependence of the gap function 
can be gauged away leaving us with an effective one¬ 
dimensional problem of the Richardson-Gaudin form [45, 
46] which can be solved exactly. 

In the absence of such a full solution, we resort 
to a time-dependent BCS variational ansatz |4/(t)) = 

n k (u k (f) + Hc(i)c k t c - k |J |0). This is equivalent to in¬ 
troducing a quasi-particle Hamiltonian 

H qp (t) = £ k 4(7 C kc ~ (7kA(t)4^4i4 + hc^j 

k a k 

( 3 ) 

and solving the equations of motion for the ex¬ 
pectation values ( 4(7 (t)cktj (t)) and (c_ k .|.(t)c k -|-(t)) 
with a time-dependent order parameter A(f) = 
J Ek7k(c -k±(t)ckf(t)) which is calculated self- 
consistently at each time and couples the different 
points in k space [47]. 

To drive the system out of equilibrium we consider 
the simplest and very popular quantum quench protocol. 
We take at t = 0 the ground state of the Hamiltonian 
[Eq. (3)] with a given value of the interaction parameter 
J = Ji and calculate the time-evolution according to 
the same Hamiltonian but with a different interaction 
J = Jf. It is convenient to discuss the results in terms 
of A j and A f , the gaps which the system would have at 
equilibrium for Ji and Jf , respectively. 

In Fig. 1 we plot the time-dependent order parame¬ 
ter A (t) for s- and d- wave symmetry and for four differ¬ 
ent values of the quench parameter Ai/Af, ranging from 




FIG. 1. (Color online) Plot of the gap dynamics for the d- wave 
(full red line) and s-wave (dashed black line) symmetries and 
for different quench parameters Ai/Af. (a) 0.001, (b) 0.2, (c) 
4.0, and (d) 5.2. Inset of panel (b): log-log plot of the local 
maxima of the gap. 


very small to very large ratios. In the s-wave case we 
recover three different dynamical regimes, in accordance 
with previous studies [38]: For small Ai/Af the order 
parameter exhibits persistent oscillations between two 
limiting values; for intermediate Ai/Af it has damped 
oscillations towards a non-equilibrium asymptotic value 
A st ^ 0; and at large Ai/Af it has an overdamped ex¬ 
ponential decay to zero. 

In the d-wave case we find damped oscillations for ev¬ 
ery Ai/Af except for very large values. On the qual¬ 
itative level, it is remarkable the disappearance of the 
regime with persistent oscillations [panel (a) of Fig. 1] 
and, more importantly, the much faster decay of the gap 
oscillations as compared to the s-wave case [panels (b) 
and (c) of Fig. 1], On the contrary, for very large Ai/Af 
the behavior is similar, with both s- and d -wave gap de¬ 
caying exponentially to zero. Indeed, as pointed out in 
Ref. [38], this regime is similar to a complete switch-off 
of the interaction, a case in which the structure factor 7 k 
has clearly little influence. 

The increased damping of the d -wave gap dynamics is 
a signature of the existence of low-energy excitations, as 
it can be understood at least in the case of small quenches 
Ai/Af ~ 1 for which we can calculate the linear response 
theory variation 6A(t) = A(t) — A* [47] 

2 2 

SA(t)/Ai oc ^ (1 - cos 2Ekit) (4) 

k 

where E ki = \/e\ + 7 k A?. The time-dependent contri¬ 
bution to Eq. (4) is dominated at long times by the low- 
energy modes, as it is evident if we replace the sum over 
momenta with an energy integral and we change vari¬ 
ables in order to introduce the superconducting density 
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FIG. 2. (Color online) Plot of the stationary gap A st /A/ 
as a function of the quench parameter A i /A f for the d- wave 
(full red line with full circles) and s-wave (dashed black line 
with empty circles) symmetries. For Ai/Af < 0.2 the s- 
wave gap exhibits undamped oscillations around the value 
indicated with a square, in this case the circles indicate the 
extrema of the oscillations. Curves without symbols: equi¬ 
librium gap A (T*) corresponding to an effective temperature 
T*. 


of states p(E) = J2 k S(E — -Ek)- For the s-wave super¬ 
conductor the density of states has a sharp edge at A, 
where it has a squared-root divergence. This leads to 
power-law damped oscillations with frequency 2 A,. The 
d -wave symmetry introduces a qualitative difference in 
the density of states, which diverges only logarithmically 
at A i and has finite value for energies down to zero. This 
results in oscillations which damp much faster. 

Steady state. - The different dynamical regimes and 
the long-time gap values for s- and d -wave symmetries 
are summarized in Fig. 2 where, following Ref. [38], we 
plot A st as a function of the quench parameter Ai/Af. 
For any value of Ai/Af we find that at long times the 
quench leads to a reduction of the gap with respect to the 
zero temperature equilibrium value A f. The difference 
between the s- and d -wave cases occurs for Ai/Af < 0.2, 
where the d -wave gap goes to a stationary value while the 
s-wave gap exhibits undamped oscillations [38]. On the 
other hand, for A*/A f > 0.2 the gap reaches essentially 
the same asymptotic value, despite the much faster decay 
of the d -wave gap. 

It is important to emphasize that in the absence of 
pair-breaking scattering terms and of any real dissipa¬ 
tion mechanism, the system persists in a nonequilibrium 
state. In particular, the expectation values (cj co .(t)ck C r(t)) 
and (c_kj,(f)ck-|-(f)) do not come to a steady state and a 
stationary value of the gap is eventually reached only as 
a result of destructive interference between different mo¬ 
menta, a phenomenon which can also be interpreted in 
terms of a quenclr-induced decoherence [48]. 

For completeness, in Fig. 2 we also plot A(T*), the gap 


for a system in equilibrium at the temperature T* corre¬ 
sponding to the energy pumped into the system through 
the quench (lines without symbols). The system could 
eventually reach this thermal value if we include scatter¬ 
ing processes not contained in the Hamiltonian [Eq. (3)]. 
While the overall behavior of A (T*) is qualitatively sim¬ 
ilar to A st , the quantitative difference is substantial con¬ 
firming the nonthermal character of the asymptotic sta¬ 
tionary state. 

Spectral features. - We now focus on the spectral prop¬ 
erties of the transient state in the d-wave case, thus mov¬ 
ing a first step towards a comparison with recent time- 
resolved ARPES experiments on high-T c cuprates [10- 
12]. These have tracked the evolution of the quasiparticle 
energy and weight Z k, which are directly related to 
the lesser Green’s function G//(t,t') = *(cLW c k <r (*'))> 
whose Fourier transform at equilibrium and zero temper¬ 
ature has a peak at negative energy 

(w) = Zk eq 5(ui + E\f) (5) 

7r H 

corresponding to the quasiparticle energy E^ with a 
weight 

^keq = (1 - g). (6) 

In principle, from the quasiparticle Hamiltonian 
[Eq. (3)] we can also calculate the out-of-equilibrium 
lesser Green’s function. However, the presence of an arbi¬ 
trary time-dependent order parameter A(f) and the fact 
that out of equilibrium the Green’s function depends on 
both time arguments, makes this a rather challenging 
task that we leave for further studies. Here, in order 
to proceed analytically and obtain some physical insight 
on the main effect of the quench, we exploit the ob¬ 
servation of the extremely fast dynamics of the d -wave 
gap. Hence, we assume that the quasiparticle modes 
do not have enough time to rearrange and approximate 
the actual dynamics with a sudden change of the gap 

A(t) = A i 0(-t) + A st 0(t). 

This approximation is also based on some experimental 
results. In particular Refs. [11, 12] have highlighted two 
clearly distinct time scales for the gap dynamics following 
the pump pulse: on a short interval of approximately 
0.3 ps the gap reaches a value smaller than at equilibrium 
(A st in our model) and on a longer time it relaxes to 
the equilibrium value, typically attained after 10 — 20 ps. 
Our sudden approximation for the gap dynamics should 
therefore be reasonable for times of a few picoseconds 
after the pump pulse. It is interesting to notice that on 
these time scales an effective temperature picture of the 
ARPES spectra is not adequate, as firmly pointed out in 
Ref. [12]. Finally, if we consider the maximum gap value 
at equilibrium for Bi2212 A f = 60meV we can estimate 
a time scale of about O.Olps for Fig. 1. 

Within this sudden approximation for the gap dynam¬ 
ics, we derive an analytical expression for G ^ (i, t') which 
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is convenient to Fourier transform with respect to the 
time difference t — t' and average over the waiting time 
t' to finally obtain [47] 

-*AG< neq (w) = Z~ neci S(u}+E k ) + Z+ ne(i d(u}-E k ), (7) 

where E k = sjz\ + q^A^. and [ are the out of equi¬ 
librium positive and negative energy weights. We no¬ 
tice that as a result of the sudden excitation also the 
positive energy peak has a finite occupation. In the fol¬ 
lowing we will focus on the negative energy peak, which 
is the one observed in ARPES, and compare its weight 
i?kncq to the equilibrium case. To this extent it is im¬ 
portant to discuss first the interpretation of the quench 
protocol in the framework of pump-probe experiments. 
In the standard picture, the quench is used to describe 
the change of a Hamiltonian parameter, in this case 
J(t) = JiO(—t) + Jf&(t), as a result of an external pertur¬ 
bation. However this scenario -which can be realized in 
cold-atom systems- is not directly relevant to solid-state 
experiments. In this context the Hamiltonian parame¬ 
ters can be considered largely independent of the excita¬ 
tion process and the quench is merely a theoretical tool 
to study the evolution of an out-of-equilibrium state. In 
this approach J t and A, are just used to parameterize the 
initial state which results from the impulsive excitation, 
whereas the interaction parameter Jf which controls the 
time evolution has to coincide with the actual interaction 
that characterizes the material. 

In this light, it is appropriate to compare the out-of- 
equilibrium spectral weight neq with the one at equi¬ 
librium [Eq. (6)] with A f. The result of this calculation 
reads [47] 

Z kneq = 1 ~ £k / E k /I sj + A st Aj7^ \ 

^keq 1 —£k/Ekf \2 2E k E k { J 

which takes a particularly clear and interesting form if 
we expand in the neighborhood of the nodal lines, i.e., 
for yk —> 0. In this case we obtain •Z' kneq /-Zkeq = 
(A st /A f ) 2 < 1 for g k > 0 while Z~ ne JZ keq = 1 - 
«7k/4fk for e k < 0, with a = A? — Af +2Ag t — 2A st Ai > 
0. In other words the nodal spectral weight is always re¬ 
duced with respect to the equilibrium value, except pos¬ 
sibly in a small region close to the Fermi surface. This 
reduction is clearly a nonthermal effect since finite tem¬ 
perature excitations of quasiparticle would lead to a re¬ 
duction proportional to the Fermi function. 

Conclusions. - We have studied the real-time dynam¬ 
ics of a simple model of d -wave superconductor excited 
by a sudden perturbation. For every value of the quench 
parameter Aj/A f the system relaxes to a non-thermal 
stationary state with a gap parameter smaller than at 
equilibrium. The presence of nodal gapless excitations 
results in a much faster dynamics compared to the uni¬ 
form s-wave case and to the disappearance of the regime 


of undamped oscillations which characterizes an isotropic 
superconductor when the quench parameter is smaller 
than 0.2. 

We have derived an analytical expression for the 
momentum-dependent photoenrission spectral weight, 
which demonstrates a strong dependence on the momen¬ 
tum. In particular we have found a strong suppression 
of the weight at the gap nodes which is clearly of non¬ 
thermal nature, a result which is consistent with recent 
time-resolved ARPES experiments. 
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PACS numbers: 

I. LINEAR RESPONSE REGIME AND ENHANCED DAMPING OF A DWAVE SUPERCONDUCTOR 

We start from the time-dependent BCS Hamiltonian 

H(t) = £k 4<r Ck CT — ^2 Ap(t) cL k ^c_p.(. c p ^ (1) 

kcr k p 


for a sudden quench of the pairing interaction Jk p (t) = ^(t)7k7p where J(t ) = Ji0(—t) + Jf0(t) and with 7 k = 
(cos k x — cos k v ) for d-wave and 7k = 1 for s-wave. The normal and anomalous Green’s functions on the Keldysh 
contour are defined as 


Gk(M') = -i{T K c\^(t)c\^{t')) 

(2) 

Gk(t,t') = —i(T)c4 k ^(i)c_k^(t)) 

(3) 

Fk(t,t') = —^(T/cCkt(^)c-k4(^ , )) 

(4) 


where Tjc is the Keldysh time-ordering, from which we can obtain the time-dependent superconducting (SC) order 
parameter Ak(t) = (c_k^(t)ck^(t)) as 


Ak {t) =-iF£(t,t) (5) 



FIG. 1: (Color online) Plot of the gap dynamics for d-wave (red, bottom panels) and s-wave (black, top panels) for different 
quench parameters Ai/Af: 0.9 (left panels) and 1.1 (right panels). Inset: plot in log-log scale of the local maxima of the gap. 
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We now compute the anomalous Green’s function within linear response in the strength of the quench 

SV(t > 0) = S J 7k7 P 4 t c - kJ .C-pt c pT (6) 

k P 

where SJ = Ji — Jf, around the initial BCS ground-state with J = Ji. We obtain to first order 


SF k (t,l/) = F k (t,lf) - Fg(t,t') = ^SJpp* / (7) 

pp' Jk 

We can evaluate the above expression within a BCS approximation, where we decouple normal and anomalous Green’s 
functions, to get 

SF k (t,t') = -(5J7 k A } ; [ dti [G k (t,ii)0 k (ti,t') + i'k(Mi)-Fk(ii,i')] ( 8 ) 

Jk 

where A; is the gap of the initial BCS state. Using Eq. (5) and summing the order parameter over k we get for 
SA(t) = ]T k 7 k <5A k (f) tw0 contributions 

5A(t)/Ai = SJ^2 7k [ dt i * [G k (t - ti)Gk(ti - t) + G k (t - h)G k (t i - t)] + (9) 

k Jo 

+SJ 7k [ dt i i [ F k(t - U)F<(U -t) + F<(t- h)F^(t i - t)] (10) 

k Jo 

The normal and anomalous retarded Green’s functions are readily evaluated in real-time 

G£(i) - -i0(t)({c kt (t),4 t (O)}} = -i0(t) (cos 2 (9 k e~ iEkt + sin 2 6* k e iEkt ) (11) 

F k (t) = -*0(i)({ckt(*),C-i4(O)}) = 9(t) sin 2# k sin E k t (12) 

while the advanced components are obtained by sending t —» — t and after complex conjugation. The lesser Green’s 
functions read 


G k (t) = i(cj ct (0)c kt (t)) = zsin 2 d k e iBk * (13) 

5k(f) = *(c-kt(0)c^ H (t)) = zcos 2 d k e lEkt (14) 

F k (t) = i(c- k 4.(0)c kt (t)) = —i sin d k cos 0 k e lEkt (15) 

where the angle 0 k is defined as usual in BCS theory 

cos 2 0 t = i(l + ^) (16) 

<17) 


with the ciuasiparticle spectrum E k = \Je\ + 7 k A 2 . Plugging the above results into the expression (9) we finally 
obtain 

2 

SA(t)/Ai = Sj\^ cos 2d k [cos 2 d k (e _l2£;kt — l) — sin 2 d k (e _l2£;k * — l)l (18) 

k 2i?k 

The imaginary part of this expression 

2 2 

Im<5A(t)/Aj = — <5 J cos2d k sin2£ k f = — Sj\^ fk 4 sin 2 E k t = 0 (19) 

V 2U k V 2U k 


vanishes identically so that the final result is purely real and reads 

2 2 

SA(t)/Ai = -^ES|(l- cos2 Ekt) 

k k 


(20) 
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In the d-wave case we can approximate the dispersion close to the nodes by its conelike structure, E k ~ yjvpk\ + v\k 2 
and using polar coordinates (E, 8 ) to express = E 2 cos 2 9 and 7k = E 2 sin 2 8/ A 2 we get 

c j />A 

SA d (t)/Ai = SA doo /Ai + 8vfvaA2 j dEE 2 cos 2 Et ~ 5A doo /A r + A d /t (21) 

In the s-wave case we can insert the SC density of states p(E ) = po8(E — Ai)E/y/E 2 — A? which has a sharp edge 
at E = Ai, to get for the approach to the steady state 

5A s (t)/Ai = 8A soo /Ai + ^p- J dE ^ E ~ E2 — cos2 Et ~ SA soo /A t + A s cos(2A t t + <j))/t a (22) 

The behavior of the time-dependent gap A (t) for a small amplitude quench, Ai/Af = 0.9,1.1 is plotted for reference 
in Figure 1. We see that for a uniform s-wave gap, small coherent oscillations persist at long times while they are 
quenched on much shorter time scales in the d-wave case. The inset in both figures show the decay of the time- 
dependent envelope of the gap function, which confirms a faster power law decay due to the momentum anisotropy 
of the d-wave gap. 


II. SUDDEN QUENCH APPROXIMATION 

In the limit of a very fast order parameter dynamics, we can assume the quasiparticle (qp) Hamiltonian to experience 
an almost sudden change of the pairing field Ak(t) = A(t)qk, with A (t) = A,;0(— t) + A st 8(t). The problem can be 
than solved analytically by a time-dependent Bogolubov transformation. We start by diagonalizing the initial qp 

Hamiltonian, that in the Nambu formalism Tk = c_k^ reads 

Hq P {Ai) = f k (£ k T* + A.;7k t x ) I'k (23) 

k<r 


in terms of a new set of fermions 


$ k = exp (i8 ki T y ) (24) 

such that 

H qp {Ai) = ^ $ k £kUA (25) 

k 

with E^i = + Ajb. Similarly for the final qp Hamiltonian, 

H qp (A st ) = ^ T k (ek t z + A st 7 k T x ) ^k (26) 

k a 


we can introduce a new set of fermions 


T k = exp (iOkTy) ^k (27) 

such that 

H qp (A at ) = Y, T k ^kuTk (28) 

k 

with Ekst = \/ £ k + Then, we can write the time evolution of the original fermionic field as 

*k(t) = M k (t)$ k (0) (29) 



hit) g k {t) \ 
Mty ) 


with 
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FIG. 2: Plot of Z keq and of the ratio Z kne(J / Z keq as a function of £k for different polar angle 9 in k-space and for two different 
quench parameter A;/A/. 


and 


f k (t) = cos E k t cos6 ki — is'mE k t cos (20k — 9 ki ) (30) 

,9k {t) = — cos E k t sin #ki — * sin E k t sin (20k — 0ki) (31) 

We can then define the following Nambu Green’s function (Gf) on the Keldysh contour, involving both normal 
(diagonal) and anomalous (off-diagonal) components 

G ka p{t, 0 = -i(Tx * ka(t)* k/»(0) (32) 

and focus in particular on the lesser Gf, G^ a g(t,t') = i{^ k p(t')^f ka (t)) whose normal component (a = f) = 1) 

G k (t,t') = *(4 t (t')c k tW) (33) 

is related to the ARPES signal, which typically probes the matrix element for transition from occupied to empty 
states. A simple algebra gives the result 

G<(t,f') = i [n(E ki ) + (1 - n{E ki ))g k (t)g* k {t')\ (34) 


where we have introduced the Fermi function n(E) = l/(e /3E + 1) with /3 = 1/T the inverse temperature of the initial 
state. It is easy to see that in equilibrium, i.e. no quench of the pairing field A,; = A/, we recover the standard result 




n(E ki ) cos 2 0 k i e lEki(4 4) + (1 - n{E ki )) sin 2 0 ki e lBki(t 4) 


(35) 


which in Fourier space gives 

G k {u) = i [n(E ki ) cos 2 0 k * S(w - E ki ) + (1 - n{E ki )) sin 2 6 kl S(co + E ki )\ (36) 

At low enough temperature, n{E k i) — 0 and only the negative energy peak is present and this will be the one whose 
evolution with the quench amplitude we will be following. In the non-equilibrium quenched case the lesser Gf will 
depend in general from both t — t! = r and t + t' . It is convenient to use r and t as independent variables and Fourier 
transform with respect to the former 

G k (w,t) = J dre lulT G^(r,2t — r) = i 

with respectively 

A k = ^n(E ki ) (cos0 k i - cos(20 k - ^k*)) 2 + ^(1 - n(E ki )) (sin0 ki + sin(20 k - S ki )) 2 


A k + °±e- 2lEkt 


6(ui + E k ) + ( B k + 


C k 


2 e 2lEkt ) 6 (oj - E k ) 


(37) 


(38) 
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B k = ^n(E ki ) (cos e ki + cos(26» k - 0 ki )) 2 + |(1 - n(E ki )) (sin6> ki - sin(26» k - 9 kl )) 2 (39) 

G k = ^n(Ski) (cos 2 6»ki - cos 2 (26» k - 0 k »)) + ^(1 - n(E ki )) (sin 2 6» kl - sin 2 (26> k - 0 ki )) (40) 

While it is tempting to identify the coefficients in front of the two delta functions as time-dependent spectral weight 
they are not positive definite unless we average over the waiting time t 

G<(w) = lim t * f dTG<{u,T) (41) 

t Jo 


to finally obtain 


(-i/n)G£{u) = Z Wneq S(u + E k ) + Z£ neq 5(uj - E k ) (42) 

where, respectively, Z^ neq = A k and Z^ ncq = B k . We can now look at the spectral weight associated to the negative 
frequency peak and compare it with the equilibrium one, for a final gap Ay. The result reads 


^kneq 
Z keq 


f 1 ~ £k/Ek \ /I ej + A^A t 7 2 \ 

\l-e k /E kf ) V2 2 E k E ki ) 


(43) 


As we discussed in the main text, we can expand the dispersion close to the nodes, for 7kA/e k < 1 to obtain an 
analytical expression 


'Z kneq 
Z keq 


(A st /Ay) 2 £ k > 0 

1 - a7k/4£k £ k < 0 


z~ 

where a = A 2 — A 2 + 2A 2 t — 2A st A, turns out to be positive for all quench parameters. In both cases we get ^ k neq < 1 
from which we conclude that, at least far from the Fermi surface the nodal spectral weight is reduced with respect to 
equilibrium. We plot the individual weights as well as their ratio (43) as a function of energy in figure 2, for different 
values of the quench parameter Aj/Ay and for different polar directions in k-space. We notice that the ratio is lesser 
than one for most values of the quench parameter, unless A, Ay and in a region close to the Fermi surface that 
shrinks approaching the nodal direction, thus confirming qualitatively the analytical result. 



